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1. Introduction 

In a previous paper (referred to hereafter as I) we have studied a symmetrical, 
semiconducting double quantum dot system, with interdot capacitive coupling in 
addition to the usual intradot Coulomb interaction. Employing the numerical 
renormalization group (NRG) approach [2 El EI an d focussing on the regime 
with two electrons on the double dot, a rich and diverse range of physical behaviour 
was shown to arise on increasing the ratio U' /U of interdot to intradot coupling 
strengths (jaj. The system first evolves continuously from a spin-Kondo state, where 
the dot spins are in essence separately quenched (SU(2) x SU(2)), to an SU(<i) Kondo 
state with entangled charge and spin degrees of freedom arising at U' /U = 1. This in 
turn precedes a rapid but smooth crossover to a charge-Kondo state in which a charge- 
pseudospin is Kondo quenched; followed by the suppression of charge-pseudospin 
tunneling, a continuous collapse of the underlying low-energy Kondo scale, and a 
Kosterlitz-Thouless quantum phase transition at a critical U' c to a degenerate charge- 
ordered state, itself a non- Fermi liquid with a In 2 residual entropy. 

Our primary emphasis in I was on the structure, stability and flows between 
the underlying RG fixed points, on the overall phase diagram and evolution of the 
characteristic low-energy Kondo scale arising for all U' < U' c in the strong coupling 
phase; and on static physical properties such as spin- and charge-susceptibilities, 
including some exact results for associated Wilson ratios. In the present paper by 
contrast we focus on dynamical properties of the system, specifically the w-dependence 
and [/'-evolution of local single-particle spectra, which at the Fermi level in particular 
determine the (T = 0) linear differential conductance across one or other dot (the 
interdot coupling is purely capacitive, so there is of course no transport through the 
pair of dots). 

A number of theoretical issues are discussed in section^ These relate in particular 
to the 'two-self-energy' description that we show is required to describe the broken 
symmetry charge-ordered (CO) phase, its connection both to the potential scattering 
inherent to the CO phase and to the 'single' self-energy of conventional field theory; 
and hence insights into the non-Fermi liquid character of the CO phase. Numerical 
results for dynamics are presented in section [21 including the important issues of 
universality and scaling of the low-energy single-particle spectra throughout the strong 
coupling phase. 

2. Theoretical issues 

We focus then on the local single-particle spectrum for dot i, Di(u) = —-\vaGu{u)), 
with Gu{lo) the Fourier transform of the (retarded) single-particle propagator Gu(t) = 
— i9(t){{ci a (t), c\ a }). In terms of the local spectrum, the zero-bias differential 
conductance across dot i at T = 0, here denoted by gt, is given by [J] 



with Ji and Vi respectively the current and voltage across the dot (and T — TrV 2 p the 
hybridization strength, as in I) . The T = linear differential conductance thus probes 
the single-particle spectrum at the Fermi level, u> = (and at finite, low bias voltage 
V s d, TrTDi(ui = eV s d) provides an approximation to the conductance). More generally, 
from the discussion given in I, one might expect the w-dependence of Di{oj) to exhibit 




(2.1) 
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Figure 1. NRG results for the zero-frequency spectral density irTDi(u) = 0), or 
equivalently the linear differential conductance gi/(2e 2 /h), as a function of U' for 
fixed U = 7. The dotted line marks the phase boundary U' c ~ 7.046. 



a rich evolution as the interdot interaction U' is progressively increased through the 
spin-Kondo regime (V < U), across the 5(7(4) point U' = U into the charge-Kondo 
regime; and then towards the quantum phase transition occurring at U' c and into the 
broken symmetry charge ordered (CO) phase. 

To motivate subsequent discussion, we shown in figure^NRG results for the linear 
differential conductance gi/(2e 2 /h) — ttTD^lu — 0) vs the interdot U' = U'/(ttT) for 
a fixed intradot U = U/(irT) = 7. As discussed further below, the transition between 
the strong coupling (SC) and CO phases (here occurring at U' c ~ 7.046) is clearly 
evident: throughout the SC phase the conductance is 'pinned' at the unitarity limit 
of <7i = 2e 2 //i, while it drops discontinuously on entering the CO phase and decreases 
montonically thereafter with increasing U'. 



2.1. Two-self-energy description 

In considering single-particle dynamics a number of rather subtle issues arise, which 
as we show have important ramifications e.g. in understanding the non-Fermi liquid 
nature of the CO phase. The first concerns the 'two-self-energy' description that is 
required to describe the CO phase. 

The SC phase is a Fermi liquid state, adiabatically connected to the non- 
interacting limit U = = V and with Kondo screening ensuring that the ground 
state is locally non-degenerate. As such the traditional 'single self-energy' description, 
based ultimately on perturbation theory in the interaction, is appropriate; where (with 
oj + = lu + i0+) Gu(uj) is expressed as 

G «H = , ^ ^7-^ (2.2) 



+ *r-Ei(w)] 



in terms of the conventional self-energy Si(w) (= Sf (u) - iS| (to) with the imaginary 
part Ef(w) > 0). 
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The CO phase by contrast is not adiabatically connected to the non-interacting 
limit (and as such is a non- Fermi liquid), as evident e.g. from the phase diagram figure 
4 of I, where the CO phase arises for U' > U' c and the critical U' c as a function of U 
satisfies U' C {U) > for all U > 0. This phase is of course characterised by charge 
symmetry-breaking, the ground state being doubly degenerate ((nL,TiR) = (2,0) and 
(0,2) occurring with equal weight in the fixed point Hamiltonian, section 4.1 of I). 
For this reason, as is obvious from its Lehmann representation, the local propagator 
in the CO phase is of form 

Gu(lj) = |[G ii;A M + G ii;B M]. (2.3) 

Here Gu ;a (ui), with corresponding spectral density Di a (uj) — —-^ImGu- t0 ,(uj), denotes 
the propagator obtained from one or other of the degenerate ground states a = A or 
B mentioned above. This in turn may be expressed as 

GuA") = ^ ~ l r ~r- ( r, (2.4) 

in terms of a self-energy £j a (u>) (= — iT,j a (ui)). The two-self-energy (TSE) 

description embodied in equations l|2.3l I2.4f> is a necessity and not a luxury when 
considering the CO phase (although equation 11'. 21 may obviously still be used to 
define a single-self-energy in the CO phase, as considered in section 2.3 below). If 
one sought for example to calculate Gu(u>) diagrammatically in the CO phase, one 
would construct self-energies from charge symmetry-broken mean-field propagators 
(unrestricted Hartree Fock); the double degeneracy of the mean- field saddle points 
ensuring that it is the Y>i a (oj) and hence Gu- a (uj) that are thereby calculated, with 
the full Gu(uj) then obtained from (|2.3|) . 

Notice also that, trivially, in a case where Hi a (u>) = is independent of a, 

the TSE description reduces formally to the conventional single self-energy description 
of equation l|2.2|l . It thus encompasses the SC phase simply as a special case. 

We now consider the implications of a TSE description, considering primarily the 
CO phase. The first point to note here is that, at the Fermi level u> = 0, the imaginary 
part of the self-energies Y<i a (u>) vanishes, 

ti a (u; = 0) = (2.5) 

for a = A,B. This reflects the fact that the fixed point Hamiltonian for the 
CO phase, section 4.1 of I, is characterised purely by potential scattering of the 
conduction electrons (a point pursued in section 2.2 below). It can also be shown 
by considering the diagrams for Y,i a (tu) constructed from charge symmetry-broken 
mean-field propagators. From equations Q2.41 12.5(1 it follows directly that the Fermi 
level spectral density is given by 

where ti a denotes the 'renormalized level' 

eto = e + E&(w = 0) (2.7) 

(satisfying the obvious 'left/right' symmetry e^ Q = e^ g with a = B or A for a = A 
or B respectively); and such that the linear differential conductance 

7rr£)i(o) = i J2 OrA«(0)). (2.8) 

a=A,B 
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This result may be recast equivalently in terms of the conduction electron phase 
shift S ia = S la (uj = 0), where S ia (uj) = arg(T ia (u;)) with T ia (uj) = U 2 G ii;ct (u) the 
conduction electron t- matrix. Using equations Q2.4I I2.5JI gives <5 JQ = arctan(^) and 

hence 7rTD ia (uj = 0) = sin (5i a ), i.e. 

7rrA(w = 0) = | sin 2 !' 5 -) (2.9) 

a=A,B 

for the linear differential conductance. 

The considerations above relating the Fermi level spectrum to the renormalized 
levels and/or phase shifts are quite general. For the particular case of particle- 
hole symmetry (e = — ^ — U' , for which the explicit NRG calculations we show 
here are performed) one has the additional, intuitively obvious reflection symmetry 
about the Fermi level, -Dm(^) = Dib{— w). For the self-energies this corresponds to 

^LlM = E iB (-<*>) and 

e + E&(w) = -[e + E? B (-w)]. (2.10) 
The renormalized levels (equation l|2.7[l ) thus satisfy 

CiA = —CiB (2-11) 

implying equal and opposite scattering phase shifts, SiA = —&iB (as noted in section 
4.1 of I), such that the linear differential conductance ttTD^Q) = sin 2 (<5j Q ). 

As already noted the SC phase is recovered simply as a limit of the above 
analysis, with E Jct (o;) = independent of a. At particle-hole symmetry in this 

case, equations 1(2.71 I2.10|) show that the renormalized level vanishes by symmetry 
(corresponding to a phase shift 6 = 5), and hence from equations 1)2. 6jl or 1(2. 10 j) that 
TrTDi(0) = 1 - the unitarity limit, seen clearly in the NRG results of figure ^ for 
U' < U in the SC phase. 



2.2. Potential scattering 

We now relate the above results for single-particle dynamics at the Fermi level to 
potential scattering. To this end consider the Hamiltonian 

H * S = E e *°L«Wa + Kia E a L<V«r ( 2 - 12 ) 
k.a k.k' ,(T 

in which potential scattering is operative in conduction channel i (with strength 
Ki a ), and for which the conduction electron t- matrix Ti a (w) is readily calculated 5 . 
Equating the resultant Tj Q (0) to Ti a (0) = V 2 Gu- a (0) considered above, one finds that 
(a) £f Q (0) = , i.e. equation 1)2.5(1 is recovered; and (b) that the renormalized level 
€i a is related to the scattering strength by 

= (2 ' 13) 

1 irpKia 

(with p as usual the conduction band density of states). 

Equation (|2.12Jl is in essence the fixed point Hamiltonian for the CO phase. More 
precisely, referring to equation (4.2) of I and noting that fl ia oc J2k a kia' the fixed 
point Hamiltonian is of form ^ i H[ s (both channels are of course involved); with 
Ki a = (n, — 1)K for the particle-hole symmetric case considered explicitly, where the 
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^A(O) = - , - „_^. 2 (2.15) 



dot charges (n/,, Tin) = (0, 2) or (2, 0) correspond respectively to the broken symmetry 
states we have denoted by a = A or B. Hence, considering explicitly the i — L channel, 

— - — - — (2 14) 

r - n P K ~ r ( ' 

with corresponding phase shifts 5iA — arctan(7rpif) = —SiB, such that from equation 
(|2.9|) the linear conductance is given by 

1 

1 + iJiTipK) 7 

in terms of the potential scattering strength pK (and is obviously independent 
of whether i = L or R is considered). It is straightforward to relate pK in 
turn to the potential scattering coupling constant K that enters the A-discretized 
NRG Hamiltonian for the CO fixed point (equation (4.3) of I): specifically, pK = 
([1 — A _1 ]/[2 In A])K. Since K itself is directly determined numerically in the NRG 
calculations, the linear conductance can thus either be determined this way from 
equation Ij2.15|l . or as the lj = limit of the single-particle spectrum Di(u>) (determined 
via the standard NRG method [H| reprised at the beginning of section 3). In practice 
we find the two methods to be in very good numerical agreement. 

As detailed in section 4.1 of I the CO phase corresponds to a line of fixed points, 
one for each U' > U' c . The critical endpoint, occurring at U' = U' c + on the CO side of 
the transition, corresponds to a critical K c given by equation (4.8) of I, and thus to 
a critical npK c = tan[§(l — -j^)]- From equation (|2.15|) the differential conductance 

gi/(2e 2 /h), which is unity throughout the SC phase for U' < U' c , thus drops abruptly 
at the transition to a value of ~ 0.197.. at U' = U' c + [Sj; as seen clearly in the NRG 
results of figure 2] This discontinuous drop in the differential conductance as the 
transition is crossed appears to be a rather general signature of a Kosterlitz-Thouless 
quantum phase transition, it being found also for a multi-level small dot close to a 
singlet-triplet degeneracy point [TO] and for a pair of Ising-coupled Kondo impurities, 
onto which maps the problem of spinless, capacitively coupled metallic islands close 
to the degeneracy point between N- and N + 1-electron states 0. 

We also add that while the existence of an abrupt drop in the differential 
conductance as the transition is crossed is not specific to the particle-hole symmetric 
case shown in figure the magnitude of the discontinuity is, reflecting the fact that 
deviation from particle-hole symmetry generates additional potential scattering of the 
same sign on the two leads. The generic form for the linear conductance follows from 
equations El EH3J> as 

7lTA(0) = -| 77^ rrr H T-r^ rrrl (2.16) 

H ' 2\l + l/(7rpK lA )2 1 + 1/(7^5)2/ 1 ' 

with the Ki a of form (again for i = L explicitly) KiA = —K + 6K and Kib = K + 5K, 
where SK represents the additional potential scattering common to each channel that 
vanishes only at particle-hole symmetry. 



2.3. Non-Fermi liquid behaviour 

The conventional description of single-particle dynamics centres on the usual single 
self-energy Sj(o;), defined by equation (|2.2(l (and the Dyson equation implicit therein). 
The analysis above has however been couched in terms of a TSE description, so the 
obvious question arises: what are the implications for Ej(u;)? 
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The general relation between S^w) and the two self-energies E iQ (w) follows 
simply from direct comparison of equations (|2.2|1 with 12. .'SI 12 .4|) , and is given by 

E i («) = 2[E ii i(w) + S iB (w)] + — — f - - • (2.17) 

uj+ + iT - [e + ±(E ij4 (a>) + 2J iB (w))J 

Since E[ q (ijj = 0) = (equation 1|2.5[) ). the imaginary part of the single self-energy at 
the Fermi level follows from equation (|2.17|) as 

E 4 '( W = 0) = [c "~^l (2.18) 

in terms of the renormalized levels Ci a , equation (|2.7I) . For the SC phase, where 
£iA — £iB generically (and = = CiB at particle- hole symmetry), equation (|2.18|l 
gives Ef (w = 0) = — just as required for a Fermi liquid state (Ef (w) oc w 2 as 
w->0). 

For the degenerate CO phase by contrast, the renormalized levels ^ ets- From 
equation (|2.18|) , the conventional self-energy thus has a non- vanishing imaginary part 
at the Fermi level lo — 0. This is of course a direct reflection of the non-Fermi liquid 
nature of the CO phase; and we emphasise that to recover this behaviour necessitates 
use of the two-self-energies Ei Q (aj) (which themselves satisfy Ef Q (a> = 0) = 0). For 
the specific case of particle- hole symmetry, = — and 



Ef( w = 0) 



' (2-19) 



{■k P K) 2 

(where equation l|2.14|l relating to the potential scattering strength is also used). 
And at U' = U' c +, using the critical irpK Cl T,j (0) is then given explicitly by 

^ {UJ T = 0) =cot 2 [f(l-^)] -4.1 :U' = U' C + (2.20) 



3. Numerical results: dynamics 

We now consider NRG results for the full frequency dependence of the local dynamics, 
analysing spectra (for a range of bare parameters) on both high and low frequency 
scales, and in particular fsection l3.3|l extracting the universal scaling behaviour in the 
SC phase. The single-particle spectrum at finite-w is obtained using the standard 
method - - a set of poles corresponding to the single-particle excitations are 
calculated from the sequence of NRG iterations, and these are subsequently broadened 
on a logarithmic scale to compensate for the initial logarithmic discretization of the 
Hamiltonian and recover the continuum. 



3.1. Spectra on high lu scales 

We begin by considering briefly the evolution of single-particle dynamics on non- 
universal 'high' energy scales — in particular the Hubbard satellites — upon 
progressively increasing the interdot interaction U' from the limit U' = where the 
dots are uncoupled. Figure |21 shows ttTD^ui) versus ui/T for a fixed U = 7 and U' = 
(dotted line), 6, 7, 8 and 9 (solid lines, top to bottom at ui/T = 10). For all U' < U, 
the Hubbard satellites are seen to be centred at approximately the same frequency, 
whereas for U' > U the positions of the satellites become [/'-dependent. 
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Figure 2. Single-particle spectra on 'high' frequency scales: 7iTDi(u;) vs u>/T 
for U' = (dotted line), 6, 7, 8 and 9 (solid lines, top to bottom) with a fixed 
= 7. 

It is straightforward to explain this behaviour by recourse to the atomic limit 
T = 0. For U 1 < U, the atomic limit ground state is (nL,n^) — (1, 1) (see equation 
(2.3) of I). This connects via single-particle excitations (for e.g. the left dot) to (0, 1) 
and (2,1) which, regardless of U', are both an energy U/2 higher than the ground 
state. The atomic limit single-particle spectra thus exhibit Hubbard satellites (poles 
in this case) at u> — ±-j, independently of U' , which is the essential origin of the be- 
haviour seen in figure|21for U' < U. For U' > U by contrast, the atomic limit ground 
state is (2, 0)/(0, 2). The corresponding excited states accessible under single-particle 
excitations, (1, 0) and (1, 2), are now an energy V — U/2 higher and thus become in- 
creasingly separated from the ground states as U' is raised; generating [/'-dependent 
Hubbard satellites at ui = ±{U' — U/2), as again seen in figure|21 

Before considering dynamics on the key low-energy scales, we make a general 
remark. NRG spectra may of course be calculated from either the even or the 
odd set of iterations (fixed points are inexorably fixed points of the square of the 
RG transformation |21I3]); and the same spectrum Di{u>) naturally results in either 
case. For the CO phase arising for U' > U' c the NRG ground state is as one would 
expect doubly degenerate, regardless of whether even or odd iterations are considered. 
The two-self-energy description considered above is thus ubiquitous in this case. For 
W < U' c in the SC phase by contrast there is an important difference between even- TV 
and odd-iV iterations, just as found in jllj for the single-impurity Anderson model. In 
the even-^V case the NRG ground state is always (for any iteration) non-degenerate, for 
which reason an effective single-self-energy description of dynamics obviously arises. 
But for the odd-TV iterations the NRG ground state is always found to be doubly 
degenerate (and with 'overlap' on the dots). This leads directly to an inherent two- 
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self-energy description of local dynamics, even in this non-symmetry-broken SC phase. 
We emphasise again that, whether even or odd iterations are considered, the resultant 
Di{uj) obtained from the NRG calculations is the same. But the key point is that a 
two-self-energy description of dynamics, which is a necessity in describing the broken 
symmetry CO phase, is equally applicable in describing the SC phase; a point not 
widely appreciated, but one which also underlies the Local Moment Approach to 
correlated electron systems ^JEIE]- The situation sketched here is directly parallel 
to that detailed in to which the reader is referred for further discussion. 

3.2. Spectra with increasing U' 

We now consider dynamics on the low-energy scales characteristic of the Kondo 
resonance that arises in the SC phase. For convenience we shall analyse the two 
regions < U' < U and U' > U separately. The former displays the crossover from 
SU(2) to SU(4) behaviour on increasing U' towards the St/ (4) point U' — U, whereas 
in the latter region we show how the Kondo resonance is destroyed as U' moves 
through the phase transition. Throughout the SC phase the Fermi level spectrum 
is (as above) pinned at the unitarity limit, -kYD^lo — 0) = 1; and we note that this 
behaviour is accurately captured by the numerics on employing the standard A^-factor 
(r — > Aj\T) [3] that corrects for the finite-A discretization inherent to the NRG. 

Figure|3fa) thus shows irYDi(uj) versus oj/Y for a fixed, strongly correlated U = 7, 
with increasing U' = 0, 6.5, 6.9 and 7 (in order of increasing resonance width). As 
expected from the thermodynamic results discussed in I, the SU(2) physics of U' = 
is found to persist with increasing U' until one gets very close to the SU(4) point. 
For U' = 6.5 for example, the functional form of the Kondo resonance is essentially 
the same as for U' — 0, the only difference being that its width is increased slightly 
in proportion to the Kondo temperature Tk (the evolution of Tk itself being shown 
explicitly in figure 5 of I) . 

As U' moves closer to U however, the LM SU< - 4 ' 1 fixed point begins to play 
an increasingly important role in shaping the low-frequency form of the spectrum. 
This is evident in the case U' — 6.9, where there is a distinct shoulder in the 
Kondo resonance when u> /Y is on the scale of the crossover from YM su ^ho IM SU ^ 
(~ (U — U')/wT = 0.1). We discuss this point in more detail when considering the 
scaling characteristics of the spectrum in section 13.31 for now we simply note that at 
frequencies below the crossover scale the spectrum is again of the same SU(2) form 
as for U 1 = 0, whereas on higher scales the tails of the resonance approach those of 
the U' = U SU (4) spectrum which is also shown in the figure. 

Figure |3Jb) illustrates the situation on the other side of the SU(i) point: we 
show TiYDi{ui) versus oj/Y, again for fixed U = 7 but now with U' — 7, 7.03, 7.044, 
7.048 and 7.1, spanning as such the transition occurring at U' c ~ 7.046. As shown in 
figure 5 of I, with increasing U' the Kondo scale Tk in the SC phase now decreases 
exponentially rapidly (following the Kosterlitz-Thouless form equation (6.3) of I), and 
ultimately vanishes at the critical U' c . The vanishing of the low-energy Kondo scale as 
the transition is approached is seen vividly in the behaviour of the Kondo resonance. 
This narrows rapidly with increasing U' (figure Ob for U' — 7.03 and 7.044), while 
maintaining the unitarity limit ttYD^uj = 0) = 1. The Kondo resonance collapses 'on 
the spot' at U 1 — U' c , and for U' > U' c in the CO phase leaves only an incoherent 
continuum around the Fermi level (with 7rrDi(0) < 0.2). 
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Figure 3. Evolution of the single particle spectrum TrTDi(uj) vs. ui/T with 
increasing U' at fixed U = 7. (a) shows the crossover from SU{2) to SU(A) 
physics, U' = 0, 6.5, 6.9 and 7 (from bottom to top); while (b) for U' = 7, 7.03, 
7.044, 7.048 and 7.1 (from top to bottom) shows the destruction of the Kondo 
resonance as the phase transition occurring at U' c ~ 7.046 is crossed. 

3.3. Universal scaling behaviour 

One of the most important results for the single-impurity Anderson model is the 
universality of the Kondo resonance, arising in the strongly correlated Kondo regime 
(U 3> 1) when the frequency axis is rescaled in terms of the Kondo scale Tk- Here, 
we examine corresponding issues of universality for the DQD model in the different 
regimes of the SC phase. 

First recall briefly the scaling behaviour of the spectra ttTD^lo) for U' = 0, 
where the single-impurity limit (at particle-hole symmetry) is recovered. For U » 1, 



Dynamics of capacitively coupled double quantum dots. 



11 




300 



Figure 4. Scaling spectra in the spin-Kondo regime. Main figure: 7iTD;(u;) vs. 
uj/T k for fixed U = 7 and U' = (top solid line), 6.95, 6.98, 6.99 (dotted lines 
from top to bottom) and 7 (bottom solid line). Inset: the collapse of the ST/ (4) 
scaling spectrum onto a common form with increasing U, ■wTDi(ui) vs. uj/Tk for 
U' = U = 7 (solid line), 9 (dashed line) and 11 (dotted line). Note that the U = 9 
and 11 spectra are virtually indistinguishable. 



spectra for different values of U collapse onto a common form (see e.g. when the 
frequency axis is rescaled as uj/Tk, with no dependence on the 'bare' parameters of 
the model. Once this 'scaling spectrum' is known, the spectrum for any particular 
U can thus be determined on an absolute frequency scale simply from a knowledge 
of Tk alone. The scaling spectrum naturally embodies the universal physics of the 
many-body Kondo resonance — the high-energy, non-universal Hubbard satellites for 
example (at \u)\ ~ ^) are obviously 'projected out' of the scaling spectrum. We show 
the U' = scaling spectrum, -kTD^uj) versus w/Tk, as the uppermost solid line in 
figure 01 and remark in passing that the analytical structure of it is both rich and 
non-trivial, see e.g. |15j . 

The U' = scaling spectrum pertains of course to the 5/7(2) Kondo model. 
One would clearly expect a similar scaling spectrum to exist at U' = U, in this case 
reflecting SU(4) Kondo physics. That this is indeed so is seen in the inset to figure 0] 
which shows irTDi(uj) versus uj/Tk for U' = U = 7, 9 and 11. From this the collapse 
onto a universal form is immediately apparent. Of course, this SU (4) scaling spectrum 
differs substantially from its SU(2) counterpart, as shown in figure 0] where the SU(4) 
spectrum is plotted as the lower solid line. 

Having shown distinct scaling behaviour at the SU(2) and SU(4) points, an 
obvious question arises: how do the scaling spectra evolve from ST/ (2) to SU(4) 
behaviour in the intervening region < U' < {/? As mentioned in the discussion of 
section rOl it is in fact the non- universal scale U — U' which controls this crossover; and 
the magnitude of U — U' in relation to Tk determines whether the scaling spectrum 
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will be essentially that of U' = 0, U' = U, or a combination of the two. For example, if 
U' is not exponentially close to U, the ratio (U — U')/Tk will be very large (recall that 
Tk is always exponentially small in strong coupling), the IjM su( - 4 > fixed point will thus 
have basically no influence on the low-frequency dynamics, and the spectrum should 
therefore collapse onto the SU(2) scaling spectrum. Similarly, if (U — U') <C Tk, then 
one would expect the scaling spectrum to be essentially unchanged from the SU(4) 
spectrum on all but possibly the lowest frequency scales (and even here the Friedel 
sum rule guarantees that nTDi(uj) — » 1 as lo/Tk — > throughout the SC phase). 
When (U — U') ~ 0{Tk) however, we do expect to see qualitatively different scaling 
spectra: here the crossover from SU(2) to SU{A) occurs on a scale of order Tk itself, 
and we would thus expect the spectra to show characteristics of both regimes. 

Our numerical results agree with the arguments given above. For U' < 6.9 
and frequencies in the range shown in figure ^ the single particle spectra plotted 
as a function of lo/Tk collapse essentially perfectly onto the U' — SU(2) scaling 
spectrum. For larger values of U' however, the ratio (U' — U)/Tk moves into the 
frequency range under consideration and hence the spectra begin to show signs of 
SU(4) tails. For example, the dotted lines in figure 01 show %TDi(u>) versus uj/Tk 
for U' = 6.95, 6.98 and 6.99, from which it is clear that with increasing U' towards 
U' = U , the SU (2) low-frequency behaviour holds over a progressively decreasing 
range of lu/Tk- For values of U' much closer still to U, the numerically obtained 
scaling spectra are virtually indistinguishable from that of the SU (4) line U' = U. 

The discussion above has focussed on the crossover from SU(2) to SU(4) 
behaviour, i.e. for U 1 < U. We now show that there is a third scaling spectrum in the 
SC phase, which describes the behaviour in the charge-Kondo regime (U <U'< U' c ) 
as the transition is approached and the Kondo scale Tk acquires its asymptotic 
Kosterlitz-Thouless form (equation (6.3) of I, vanishing exponentially rapidly as 
U' — > U' c ~). The main part of figure [5] shows TrTDi(ui) versus lo/Tk for fixed U = 7 
and U' = 7.03, 7.04 and 7.044: as U' approaches U' c (~ 7.046), the Kondo temperature 
Tk rapidly diminishes and the rescaled spectra collapse onto a common form. The 
scaling spectrum shows a distinct 'charge-Kondo' resonance, the tails of which tend 
not to zero as lo/Tk — > oo, but to the value sin 2 [ : |(l — -^=)] ~ 0.2 which is the 
zero- frequency value of the spectrum on entering the CO phase (section |2.2| I. 

In the inset of figure [SI we show further that this scaling spectrum is independent 
of the value of U (for strong-coupling U 1), the same spectrum being obtained close 
to the phase boundary for U = 7, 8 and 9. For each U, the value of U' has been chosen 
sufficiently close to the phase boundary that the Kondo scale Tk is well described by 
its asymptotic Kosterlitz-Thouless form (see e.g. figure 5 of I); the U's being chosen 
in practice such that (U' - U)/(U' C -U) = 0.86. 

To summarise, we have shown that the scaling behaviour of the spectra 
throughout the SC phase can for the most part be understood quite simply in terms 
of the three scaling spectra obtained in the limits U' = 0, U' = U and U' = U' c —. 
The U' — scaling spectrum in fact captures the behaviour in the vast majority of 
the phase, being the appropriate description for all U — U' ^> Tk- As one moves 
closer to the SU(A) line, the spectra at low-frequencies cross over to the SU (4) scaling 
spectrum shown in figure 01 and then beyond U' = U the spectra rapidly collapses 
onto the charge-Kondo scaling spectrum shown in figure |SJ 
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Figure 5. Scaling spectra in the chargc-Kondo regime. Main figure: 7iTD;(a;) 
vs. u)/T K for U = 7 and U' = 7.03 (dotted line), 7.04 (dashed line) and 7.044 
(solid line). Inset: ■nYDilu)) vs. ui/Tj^ for U = 7 (solid line), 8 (dashed line) and 
9 (dotted line) where ({/' - U)/(U' C - U) is fixed at ~ 0.86. Note that the spectra 
for the largest two values of U are essentially indistinguishable. 

4. Discussion 

Finally, we comment briefly on symmetry breaking perturbations which destroy the 
stability of the CO fixed point. Particle-hole asymmetry does not of course fall into 
this category. As we have emphasised several times here and in I the particle-hole 
symmetric case (e = — ^ — U') is entirely representative, and the essential physics 
robust to departure from it. However the model we have considered is 'left/right' 
symmetric: equivalent dots (e.g. cr = = e), and equivalent coupling to the 
conduction channels/leads (T^ = Tr = T). If L/R symmetry is broken it is readily 
shown that while the SC fixed point remains stable, additional relevant perturbations 
arise for the CO fixed point and hence render it unstable. Flows in the vicinity of the 
CO fixed point thus ultimately flow away from it under renormalization, and cross 
over instead to the SC fixed point. 

Since the quantum phase transition between the SC and CO phases is strictly 
destroyed by breaking L/R symmetry, the question arises: do visible remnants of the 
transition nonetheless remain? Indeed they do. To illustrate this we consider explicitly 
the case of detuning the dot levels e.g. by applying a different gate voltage to each 
dot, taking e L ^ R = e±^Se (with e = — (^ + U') as usual). Here one expects intuitively 
that the crossover scale for flows away from the CO fixed point will be determined by 
5e; and that if Se is not large compared to the SU(4:) scale T^- U ^ 4 \ then the erstwhile 
transition will be 'smeared out' but still in essence visible as a clear crossover. That 
this is so is illustrated in figure 6, where we show NRG results for the differential 
conductance 7rri?L(0) vs U' for three different 5c, viz Se/T^- U ^ = 0,0.3 and 2 (with 
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Figure 6. Linear differential conductance nTDr,^ = 0) vs U' for a fixed U = 7 
and Se/T^y^ = (dotted line), 0.3 (dashed) and 2 (solid). Note that the lowest 
U' shown in the figure is U' = 6.9. 



Tk U(4) /T ~ 0.01 the U 1 = U = 7 SU(4) scale for 5e = 0). The pristine transition 
occurring for Se = is as anticipated smeared out, on a scale determined by 5e; albeit 
remaining clearly visible (note that the lowest U' shown in figure 6 is U' — 6.9), and 
as such robust in an obvious sense. 
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